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Summary findings 

Rural roads "promote economic development but also 
facilitate deforestation. To explore the tradeoffs between 
development and environmental damage posed by roaH 
building, Chomitz and Gray develop and estimate a 
spatially explicit model of land use. This model takes 
into account location and land characteristics and 
predicts land use at each point on the landscape. 
They find that: 

■ Market access and distance to roads strongly affect 
the probability of agricultural use, especially for 
commercial agriculture. 

• High slopes, poor drainage, and low soil fertility 
discourage both commercial and semi-subsistence 
agriculture. 

• Semi-subsistence agriculture is especially sensitive to 
soil acidity and lack of nitrogen (confirming 
anthropological findings that subsistence farmers are 
shrewd judges of soil). 

Spatially explicit models are analytically powerful 
because they exploit rich spatial variation in causal 
variables, including the precise siting of roads. They are 
useful for policy because they can pinpoint threats to 
particular critical habitats and watersheds. 

This model is a descendant of the venerable von 
Thunen model. It assumes that land will tend to be 
devoted to its highest-value use, taking into account 
tenure and other constraints. The value of a plot for a 



panicular use depends on the land's physical productivity 
for that use and the farmgate prices of relevant inputs 
and outputs. A reduced-form, multinomial logit 
specification of this model calculates implicit values of 
land in alternative uses as a function of land location and 
characteristics. The resulting equations can then be used 
for prediction or analysis. 

The model was applied to cross-sectional data for 
1989-92 for Belize, a forested country currently 
experiencing rapid expansion of both subsistence and 
commercial agriculture. A geographic information system 
was used to manage the spatial data and extract variables 
based on a three kilometer sample grid. 

Three land uses were distinguished: ''natural" 
iKgetation, comprising forests, woodlands, wetlands, and 
savanna; semisnhsistence agriculture, rnmprising 
•raditional milpu (slash-and-bum) cultivation and other 
nonmechanized cultivation of annual crops; and 
commercial agriculture, consisting mainly of sugarcane, 
pasture, citrus, and mechanized production of com and 
kidney beans. 

Two dimensions of distance to market were 
distinguished: the distance from each sample point to the 
road, and on-road travel t<me to the nearest town. Data 
on a wide variety of land and soil characteristics were 
also used. 
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1.0 Introduction 

The loss of tropical forests is a major environmental concern, because it threatens biodiversity, 
contributes to global wanning, and has local effects on erosion, flooding and possibly climate. 
Considerable attention has been devoted to deforestation, and progress has been made in 
understanding its causes. But the relative impact of different causes, and the manner in which 
they interact, are poorly understood. In general, our understanding of the quantitative 
dimensions of deforestation processes remains extremely weak. 

This is especially true for a key policy issue: road-building. Roads are closely associated with 
deforestation in many parts of the tropical world. However, construction of rural roads has 
traditionally been one of the most important tools for economic development, a tool which 
moreover tends to favor the rural poor. It is Ih^erefore important to quantify the impact of road- 
building on both deforestation and development in order to assess the severity of the trade-ofif 
between environmental pr^rvation and economic growth. 

Do all forest roads cause equally severe deforestation? Whom do these roads benefit? A priori, 
we might cxjpcct deforestation to be less where population densities are low, where soils are 
unsuitable for cultivation, and where markets are distant If indeed road impacts are modulated 
by local conditions, then it may be possible to site roads so as to spur development ^ndiile 
minimiang induced deforestation. In particular there may be strong implications for choices 
between extensification vs. intensification of the road network. 

In order to explore the issue of road impacts, this paper develops and estimates a spatially 
explicit mode! of land use — one tiiat takes location and land characteristics into account, and 
predicts land use at each point on Ifae landscape. Spatially e}q>licit models are apt for two 
reasons. First, they exploit rich spatial variation in causal variables — variation wiiich is 
obscured in aggregative data (e.g. district-level means). Second, location matters. In general, we 
are interested not just in the areal extent of deforestation, but the degree to which it affects 
oitical habitats and watersheds. 
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The model presented here utilizes spatially disaggregate data, controls for a wide variety of land 
and soil characteristics, employs multiple land use categories, and is embedded in an economic 
framework. The framework, derived from von ThOnen's famous model, can be extended to 
explore the impacts of changes in commodity prices. 

The model is estimated on data describing Belize in 1989. Still mostly forested, Belize is of 
great conservation interest because of its rich biodiversity and because of its relatively large 
tracts of contiguous forest. Despite its small size, Belize exhibits a diverse array of deforestation 
processes, including encroachment by swidden agriculturalists, and forest conversion to pasture, 
sugar, citrus, and large mechanized &nns. Belize also enjoys superb documentation of land use 
and land characteristics, fecilitating Has kind of study. 

The plan of the paper is as follows. We first summarize the principal policy issues and existing 
analyses. A nontechnical presentation of our land use model follows. OBconometrically oriented 
readers are referred to Appendices A and B.) Next, we discuss tiie relevance of Belize to the 
issues, and review the main features of land use there. The results of model estimation are then 
presented, and implications discussed. The paper concludes with a summary and discussion of 
next steps. 

2.0 Deforestatton, roads, and land use: issues and anafyses 

The principal purpose of this paper is to examine the impact of roads on deforestation and land 
use. Qualitatively, the impacts are clear. New roads offer market access for timber and for 
agricultural products from previously remote areas. Roads also decrease the cost of migration, 
access, and land clearing for subsistence &rmers. In sum, road construction into forested areas 
unambiguously increases the incentives to log those territories or convert them to other uses. 
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But much hinges on the quantitative magnitude of those incentives, which we hypothesize to 
vaiy systematically over the landscape. Consider the following issues in regional and 
environmental planning: 

Rotul extensification vs. intensification: Schneider (1994) and others have suggested that road- 
building should be intensive rather than extensive. That is, road development should stress the 
creation of dense road networks around market centers rather than the extension of roads into 
low-population density areas with good soils. While this seems to be a reasonable proposition, 
we lack information about the relative oivironmental costs and development boiefits of the two 
strategies. Is road intensificaH jn a "win-win" strategy — that is, does it boost output and reduce 
environmental damage compared to an extensification strategy? What are its distributional 
implications? . The significance of these questions is underlined by the rapid expansion of the 
road network in the tropical world. Over the 1980's, Brazil's paved road network grew fiom 
87000 to 161,500 kilometers, and Indonesia's fiom 56,500 to 1 16,500 kilometers (World Bank, 
1994). In sub-Saharan Afiica, a recent review concluded tiiat "the present rural road network ... 
needs to be increased xsp to tenfold if the Ml agricultural potential of the region is to be realized." 
(Riversonero/., 1991). 

Environmental impacts assessment of forest and mining concessions. According to some 
(Kummer and Turner, 1994; Bruce and Cabarle, 1993), logging's indirect impact on deforestation 
may be greater than the direct impacts of timber removal and collateral damage to standing stock. 
Further damage results as logging roads and operations £u;ilitale access by follow-on settlers, 
who convert the logged-over forest to pasture, ])ennanent crops, or shifting cultivation. Hence 
plans for the sustainable management of forest concessions need to go beyond purely 
silvicultural considerations. Predictive models of follow-on settlement could be en^loyed in 
environmental impact assessments of proposed logging concessions, and of mining concessions 
which entail road-building. 

Conservation planning: It is expensive to set up and maintain protected areas. Conservation 
plannos have long recognized the need for an index of the threat of convosion as an aid to 
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prioritizing candidate areas for protection. Such indices have been constructed as ad hoc 
functions of population density, prior conversion, and so on. Without beliavioral grounding, 
these indices may however not be vciy accurate. Past deforestation rates, for instance, may be 
poor predictors of current rates if prices or policies luive changed. Cross-sectional variation in 
population densities usually reflects differences in soil quality, and may bear no relation to 
incentives for deforestation. Much more desirable would be a methodology which measures 
deforestation incentives in economic terms. This could be used to calculate the opportunity cost 
of conservation progruns, and could inform the design of conservation policies which seek to 
alter these incentives. 

In sum, quantitative models of land use and land use change could be used for a variety of 
envirormiental plaiming purposes. But to be useful, quantitative studies must meet several 
criteria. Firpt, they must use spatially disaggregate data. They must incorporate a wide range of 
land-use determinants, but recognize that population distribution, road placement, and land use 
change are jointly determined. Finally, they must be based on an economic framework. 

No existmg study for developing countries meets all of these criteria' Liu et al (1993) present a 
spatially detailed description of deforestation in the Philippines. They partition the landscape 
into 17 classes, based on distance from the existing road network in 1941, and show that there is 
a strong inverse relation between distance and proportion of forest lost over the period 1934- 
1988, up to a distance of 16.5 miles; from 16.S to 25 miles, the deforestation rate was constant 
This si^gestive bivariale relationship, however, may arguably reflect the action of other 
variables. For instance, if roads were first built in areas with the most fertile soils, or if road 
density were higher near emerging cities, then tiie results may reflect the effect of soil quality and 
market proximity rather than the true impact of roads. 

Reis and Margulis (1991) estimate a multivariate regression model of deforestation using county- 
level cross-sectional data from the Brazilian Amazon. The log ratio of deforested to forested 



1 For a conqnehensive survey of defoFestadoa models, see Lambin (1994). 
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areas is regressed on population density, cattle density, crop area, logging activity, road density, 
and distance to the state capital. Road density is found to be highly significant; the coefficient on 
distance to the state capital is not significantly different firom zero. Tlie equation is hard to 
interpret, though, because most of the explanatory variables are in fact jointly determined with 
deforestation. That is, we would expect land conversion processes to simultaneously result in 
deforestation, intensification of the road network, and expansion of cropland and pasture. 

An alternative modeling strategy emphasizes prices as the fundamental drivers oi° land use 
change. Multivariate regression models in this vein include Panayotou and Sungsuwan (1989) 
for forest area in Thailand and Barbier et al. (1993) for cultivated area in Mexico. Both studies 
use state-level panel data. Explanatoiy data in both include agricultural prices, provincial 
income or income per capita, population or population density, road density. Panayotou and 
Sungsuwan find that the elasticity of forest cover with respect to road density is just -O.Il . Since 
roads affect prices, however, we would expect that some road impacts are captured in the 
agricultural and timber price variables, which are found to be important. By construction, 
however, this model will attribute to prices any road impacts vdiich operate via changes in 
agricultural or timber prices; the latter is found to be quite important The Thailand study also 
finds that the elasticity of ibrest cover with respect to distance firom Bangkok is .70. Barbier et 
al., in contrast, find that increased road densi^ is associated with decreased cultivation area (they 
are not able to directly measure forest cover). Here road density may be capturing some aspects 
of urbanization, rather than changes in farmgate prices - an example underlining the need for 
higihly disaggregate geographic information. Once again, virtually all the explanatory variables 
in both studies are arguably endogenous. 

« 

Southgate, Sierra, and Brovra (1991) examine deforestation in Ecuador with a reduced-form 
model similar to that presented in this paper, but use canton-level data They regress ^ricultuial 
population on urban population, road lengtii, and soils. In a separate equation, deforestation is 
regressed on agricultural population and tenure security. Deforestation is strongly linked to 
agricultural population, which in turn is positively correlated with soil quality and witii road 
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length. The model does not however control for the endogeneity of agricultural population, or 
allow for the influence of roads on deforestation independently of their effect on population. 

3.0 A spatial model of land use 

Deforestation is just one aspect of a general model of land use. The model to be estimated traces 
its ancestry back to von Thtlnen (1826). Theoretical variants have been presented by von 
Amsberg (1994), Schneider (1994) and Hyde et al (1993). The only econometric 
implementation of which we are aware is a US application by Alig (1986). Fox et al (1994) 
present a formally similar econometric model of crop choice by Thai farmers, but without appeal 
to an economic moHel. 

The basic idea is straightforward. (See Figure 1). There is a potential rent (farmgate value of 
ou^ut minus costs of inputs) attached to each possible use of each possible plot of land. The 
model predicts, simply, that land will tend to be devoted to the activity yielding the highest rent 
Thus, in the classic example, farmers near a city find vegetables more profitable Han grain. But 
because of perishability, it is more expensive to transport vegetables than grain. 



HGURE 1 - TEIE MODEL: VON TEIUNEN REVISITED 



RENT 
$ 



Vegetables 



(expensive to 
tiSBosport) 




Land Use 



Vegetable 



Grain 



Undisturbed 
Forest 



DISTANCE TO MARKET 



A Spatial Model of Larui Use in Belixe 



7 



At some distance from the city grain becomes more profitable thon vegetables. Therefore, at a 
greater distance, where it is not economically feasible to transport any crop to the city, the land 
may be used by subsistence fanners. Further yet, the land may be undisturbed under it^i original 
forut cover. 

The gist of the model is as follov/s. (A fonnal derivation of the model and its econometric 
implementation are given in Appendix A). The potential rent associated with devoting plot i to 
u ie or commodity k is' : 

(1) Rent fork at i" (Farmgate price of k at i)*(Per-hectare production of k at i) - 

(Farmgate price of inputs to k at i)*(Qwmtity of inputs to k at i) 

Unfortunately we observe none of these variables. We do however observe the determinants of 
price and of productivity and can therefore formulate a reduced form model'' . First, following 
von ThCtaien, we assume that farmgate price of outputs decreases with distance fiom mariceL We 
assume that input prices increase with distance fiom market This should certainly be true for 
manu&ctured inputs such as fertilizer. Labor costs, too, should be relatively higher at greater 
distances fiom maik<^ coiteis, where population density is lower and woricers require 
compensating differentials to reflect the paucity of amenities (healdi, education, entertainment).^ 
We distio^sh two dimeissions of distance: the distance to the nearest road, and the travel time 
along that road to the nearest market For each commodity we can thai specify two fimctions: 

(2) Fanngate price of k at i " (distance to road at i, on-road travel time to market) 
Farmgate price <tf inputs to k at i * C|((dislance toroad at i, (m-road travel time to market) 



2 This assumes a static firamework. Dynamic issues are discussed below. 

3 Note tbat all die right hand side variables are jointly deteimiiied widi land use in a qntial equdlibrhun and 
dierefore endogenous to die model. Hence, to estiinaie(l) directly, we would faavs to instruinentiliese variables 
via auxiliary equations paraOd to diose used to derive die reduced fonn. 

4 On die odierhaiid. higher iniian costs of living ooiddresuU in a negative gcadiem of W4ges as di 
town increases. 
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The productivity of plot i for commodity k will depend on land characteristics at i. These 
include the nutrient content of the soil, acidity, workability, drainage, susceptibility to flood, and 
slope: 

(3) Per-hectare productivity of k at i « S|c(land characteristics at i) 

With suitable assumptions about functional form (see Appendix A), equations (2) aiid (3) are 
substituted into (1), yielding: 

(4) In Rent for k at i oo)^ "^^IV. distance to road + a 2k on-road time to market + a 31^ (measure of 
agricultural suitability fork) + ...+ uj|( 

Given the assumptions, the coefficients on distance are unambiguously predicted to be negative, 
and those on land characteristics will be consistent with the characteristic's predicted effect on 
productivity. Note that a random disturbance term, u, has been added. This represents tiie efifect 
of unmeasured variables. 

Finally, we add the assumption: 

(5) Plot i is used for commodities k if: rent for k at i> rent for any otfao- commodity at i 
This is a strong assiunption, witii the following underpinnings: 

Lcmd use is reversible. If cultivated land becomes uneconomic (due to a drop in crop prices, 
say), it will revert to natural vegetation. This is a defensible assertion, even in the short nm, as 
long as "natural vegetation" is broadly defined. If "forest" were distinguished as a separate land 
use, reversibility certainly would not hold in the short run, dnce abandoned land takes years to 
return to forest For many but not all forest areas, reversibility night be a reasonable assumption 
in the context of a long-nm, static equilibrium model. Large portions of Belize's forest, for 
instance, were in the recent past levelled by hurricanes. Studies in the Amazon, too, have shown 
that most abandoned plots quiddy revert to forest ^otan et aL 1994), though this is less true of 
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plots -vAach have been intensively used or scraped by bulldozers (Nepstad et al. 1991). It is 
important to stress, though, that deforestation is only one fonn of forest degradation. Regrowth 
of forest cover, for instance, does not necessarily imply maintenance of original biodiversity 
levels. 

Tenure as a determinant of rent. Equation (5) assumes that landowors will either adopt the 
highest-rent land use, or rent or sell the land to someone else who will do so. But, as Schneider 
(1994) and Hyde et al. (1993) have stressed, returns to different land uses depend strongly on 
tenure. On the fix>ntier, where land rights are poorly defined and difficult to defend, it may not 
be profitable to invest in perennial crops. Given tenure security, however, perennials may 
represent the highest value use of the land. Similarly, largeholders may refiain fi'om renting out 
land to sharecroppers, even where the latter enjoy higher returns, if land tenure might thereby be 
jeopardized. Hence it is desirable to use the land's tenure status as an explanatory variable in 
equation (3). 

Future price changes are controlled for. If today's land use affects tomorrow's land productivity 
or tenure, then the fiiture price-paths of alternative products will affect current land use decisions. 
A classic example is deforestation as a means of asserting land rights in an area where land prices 
are expected to rise. (Schneider 1994) Ideally, then, the expected path of future prices should be 
included as explanatory variables.Given (5), and assumptions about the statistical properties of u, 
statistical methods can be used to find the coefficients for (4) which best explain observed 
patterns of land use. One convraiient, but very strong, set of assumptions (see Appendix A) leads 
to a multinomial logit model, which here yields a very simple formula' : 

(6) Predicted probability that plot i is used for 

k = (predicted rent for k at i)/( Sum of predicted rents for all possible uses j at i) 



S For tedmical reasons, tbe coefficients of (4) for one land use aie normalized to zero. 
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Equivalently, under lids fonnulation equation (4) describes not just the In rent for k, but also the 
hi odds that plot i is devoted to k relative to a baseline use (e.g., undisturbed vegetation.) 
The multinomial logit model, however, requires that the unobserved effects on the rent for 
commodity k be independent of the unobserved effects for other commodities at the same point. 
This is implausible; immeasuied aspects of soil fertility, for instance, may have similar effects on 
a variety of alternative crops. In fiiture work, a multinomial probit formulation will be £^plied. 
This allows for correlation among the unmeasured effects, but is computationally much more 
conq>lex and can handle only a limited number of alternative uses. 

3.1 Dealing widi endogeneity 

We would like to interpret tiie results of tiie model as telling us the effect of tiie road network on 
agricultural land use. This interpretation would be straightforward for roads whose placement 
was not motivated by agricultural development proq>ects. For instance, some roads are installed 
for political reasons, or to provide access to a mine site, or to connect distant cities. 

In general, however, road construction and routing may be endogenous — tiiat is, influenced by 
agricultural development considerations. If loads tend to be preferentially routed tiuough 
agriculturally suitable areas, and if some aspects of suitsibiljiy are not observed, then the model 
will tend to overestimate the effect of distance fiom tiie road. A plot of land may be 
undeveloped not because it is &r from the road; it may be far fiom the road because it is not 
suitable for development 

There are two solutions to this problem. One, used here, is to employ data vtdiich provides a rich 
set of variables to control for agricultural suitability. By explidtiy controlling for the most 
important determinants of agricultural suitability, potentisd bias is greatly reduced. 

Ideally, it would be desirable to instrument the distance variables. In Mas api»oach, distance to 
road, and on-road distance to market are modelled as fiinctions of instrumental variables ^j^ch 
are maintained to be indepradent of unmeasured aspects of agricultural suitability. Predicted 
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values of the distance variables are then used in equation (4) to predict land use. Because the 
predicted values are, by hypothesis and construction, "purged" of endogeneity, the resultant 
estimates are unbiased. This strategy is appealing but difficult to implement because of the 
difficulty of Ending appmpriate instrumental veuriables. Appendix B discusses the 
instrumentation strategy fiirther, reviewing one ineffective approach and nominating a promising 
alternative for future experimentation. 

3.2 Spatial autocorrelation 

Land characteristics are quite likely to be characterized by spatial autocorrelation: places that are 
close by will tend to have similar soil types, rainfall, and so on. Some of these chaiacteristics 
may not be observable. The resulting spatial autocorrelation of disturbance terms (the in 
equation 4) is the two-dimensional analog of more ^miliar patterns of autocorrelation in time- 
series models. In general, it results in inefficient parameter estimates and inaccurate measures of 
statistical significance. 

Methods exist for dealing with spatial correlation in linear models (Anselin 198S). A 
computationally complex algorithm has been proposed for binary logit models (McMillen 1992) 
but no equivalent exists for the multinomial logit or probit models considered here. If however, 
we assume that the combined effect of the unobserved variables varies smoothly over the 
landscape, tiien a spatial trend can be used to remove autocorrelation. The two-dimensional 
equivalent of a time trend, this consists simply of a polynomial in latitude and longitude. 

4.0 Land use in Belize: context and relevance 

Belize is a small country, with about 200,000 inhabitants and about 22,000 square kilometers of 
land area. As Map 1 (IBRD 17093111) shows, only about 12% of tiie land area has been 
converted to agriculture or setflements; 65% is imder broadleaf forest, the remainder consisting 
mostiy of swamp, pine forest, and mangrove forest. (See Map 2; Table 1 gives a more 
disaggregate breakdown of land cover) 



IBRD 17093RI 
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Table 1: Land use in Belize, 1989/92 



Percentage of 


Land Use 


Area (Ha) 


Total Land Area 


Uiban/industrial 


8355 


0.38 


Sugar-cane 


64072 


2.95 


Cftnis 


12963 


0.60 


Annual crops (medianiied) 


40246 


1.85 


Bananas 


2058 


0.09 


Mango 


16S4 


0.08 


Cooca 


202 


0.01 


Cashew 


34 


0.00 


Pasture 


3S622 


1.64 


Other commercia] crops 


79815 


3.67 


Milpa fanning 


37162 


1.71 


Annual crops (non-medi.) 


19619 


0.90 


Semi-subsistence 


56781 


2.61 


Broadleaf forest (inc 


1429921 


65.77 


regrowth) 






Pine forest 


64881 


2.98 


Forest 0nc. r^nnvfli) 


1^04802 


68.75 


Savannah 


192192 


8.84 


Thicket 


84773 


3.90 


Herb, and shrub 


18844 


0.87 


Bamboo/riparian 


11518 


0.53 


Coastal strand vegetatKMi 


Z4al 


nil 

U. 11 




31255 


1.44 


Saline swamp 


34460 


1.58 


Marsh swamp 


41931 


1.93 


Bare ground 


765 


0.04 


Savannah, thicket, wetland 


41S218 


19.24 


Water bodies 


39181 


1.80 


Total 


2174187 


100.00 



Despite its small size, Belize is a reasonable exemplar of many of the issues and ciicumstances 
of deforestation. First, its population/forest area ratio, while low, is of tiie same ordo: of 
ni£^nitiide as a number of important forest regions, including Bolivia, the Congo, Gabon, Papua 
New Guinea, and the states of Amazonas and Para (See Table 2). 
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Table 2: Population compared to forest 


cover 








C.1990 


1990 


RegjoB 


People/ 


Population 


Forest area 




km2 forest 


thousands 


000 ha 


Amazonas 


1.3 


2089 


154842 


Belize 


9.5 


190 


1996 


Bolivia 


14.5 


7170 


49317 


Brazil 


zo.o 


149040 


561107 


Congo 


11.2 


2230 


19865 


Gabon 


6.4 


1160 


18235 


Papua New Guinea 


10.8 


3880 


36000 


Pari 


4.6 


S08S 


110073 


Suriname 


2.8 


420 


14768 



Sources: 

WBI, World Resources 1994-95 

except: Brazilian state populations (1991): Statesman's Yearbook 1994/95 
Brazilian state forest cover: FAO FORIS database 



Second, Belize is highly biodiverse, boasting for instance 528 bird species, as compared to 650 
for the U.S. (WRI 1994). While endemicity (the number of unique species) is low, Belize's long- 
run conservation importance is linked to its large tracts of contiguous forest (important for wide- 
ranging species such as the jaguar) and to the rapid rales of deforestation in neighboring 
countries. Third, the forests of Belize are facing increasing pressure fix)m agricultural 
development The past twenty years have seen a significant influx of refugees (up to 25,000 
according to King et al., 1993) from £1 Salvador and Guatemala. These immigrants, viho 
primarily practice subsistence agriculture, have settled around the Maya mountain block, 
significantly increasing pressure on land and forest resources in these areas. Recent immigration 
across the Westem border from the Peten region threatens the integrity of the Chiquibul National 
Park and the adjoining Forest reserves. At the same time, a booming citrus market and increased 
road construction have encouraged the spread of commercial agriculture away fit>m the 
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processing facilities. Fourth, because much agriculture is directly or indirectly supported, policy 
changes could increase or relieve pressures for agricultural conversion. 

To set the modeling work in context, we briefly review below the main features and trends in 
land use and land tenure. 

4.1 Agricultural land use^ 

Sugar. The largest single agricultural land use is sugar, at 51,000 ha ' Unusually for the 
Caribbean, it is grown mostly by smallholders. Sugar has been cultivated in northern Belize since 
the mid ISSO's, when refugees from the Yucatec Caste Wars in Mexico initially provided die 
labor supply which made cultivation possible. Cultivated area e}q>anded rapidly wilfa the 
establishment of a large processing plant in the 1950's and a second in 1967. Belize is a high- 
cost sugar supplier. Per hectare yields are low, and the product must be barged along die coast to 
Belize City, where it is transferred to ships berthed of&hore. ^Belize City lacks a deep-water 
harbor). The industry is supported by a variety of preferential trade agreements. 

Pasture. Approximately 33000 ha are devoted to pasture, by both small and large holders. Most 
beef is for domestic consumption. The industry enjoys tariff protection. 

Mechanized farming . Largeholders, many of whom are Mennonite immigrants, devote about 
41000 ha to mechanized annual crops, primarily maize, rice, and red kidney beans. 

Miilpa. Milpa is the local form of slash-and-bum agriculture, practiced primarily by ^4^yan 
Indians and by Hispanic immigrants. About 40,000 ha are under active milpa cultivation, 
primarily in the southern part of the country. The total area affected is much larger, however. 



6 Tbis section draws heavily on King et aL (1993). 

7 Land use figures in this secdon are from the 1989/92 land use survey and may differ from Minktiy of 
Agriculture rqmrts. 
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because the land is rotated under forest fallow. Maize is the most imporlant crop, produced 
primarily for self-consumption. Rice is also important Although the price is supported by the 
Belize Marketing Board, fmancial analysis suggests that the imputed wage rate for dryland rice 
cultivation is below the going W£^e for unskilled labor (even assuming zero returns to land) 
(King et al., 1993, p. 130). Milpa cultivation is likely to increase rapidly over the coming years if 
high rates of inunigration continue. This could result in increased competition for land among 
subsistence farmers, and encroachment into protected areas. 

Other non-mechanized annual crops. These cover £^proxiniately 20,000 ha and are primarily 
located near towns. It is possible that there is cross-classification error between these orops and 
milpa cultivation; in some cases a l^rmer may cultivate both kinds of crops. 

Citrus. The 1989/92 land use survey identified about 13,000 ha under citrus: too small an area to 
be accorded a distinct class in the land use model estimated below. However, citrus cultivation is 
growing rapidly. The area under cultivation is estimated to have more than doubled between 
1985 and 1990 (King et cH., 1993, p. 33). Cultivation is expected to increase further vnth the 
upgrading of the Hummingbird Highway ^elmopan-Middlesex) and the Southern Highway 
(Dangriga-Punta Gorda). While Belize has enjoyed preferential access for citrus concentrates to 
the US market, its competitive position may be threatened by Mexico under NAFTA. Lack of a 
deep water harbor and a lower scale of production place Belize at a disadvantage. 

Other export crops. Otiher crops, vAnle of potential economic importance, occupy comparatively 
littie area. Bananas occupied about 2300 ha; mango, cocoa, and cashew another 1900 ha. 

Forestry. Belize has exported timber for more tiian two centuries. Interestingly, thougih, logging 
has had littie significant impact on forest area. This is because loggmg has always focused on 
selective extraction of very higjli value species: logwood (used for dyes), mahogany and cedar. 
Extraction has tended to employ relatively low-impact methods. Railways and roads built for 
timber extraction did not caiaiyze agricultural development; in many cases roads and rsul lines 
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have reverted to forest Logging activities have, however, severely depleted the stock of 
mahogany and cedar (Alder, 1993). 

4J1 Land tenure 

Land trauire in Belize has a strong effect on land use. Table 3 shows a breakdown of Belize's 
land area according to tenure. About 22% of the country is National land, state land inherited hy 
die government in the late nineteenth and early twentieth centuries following the Mlure of die 
large private logging estates. This land is available for lease and eventual purchase and has 
traditionally been viewed as a land bank by successive governments. This is being gradually 
disbursed in line with constitutional guarantees concerning die rights of Belizeans to both own 
and 'enjoy fully* a piece of land. Unleased portions of notional land are thought to be more 
susceptible to extralegal conversion (for instance to subsistence agriculture) than other areas. 
Another 40% of the land area is privately owned land and is less likely to cxpeneace squattii^. 
About 14 % of Belize is currently protected as a national pack, wildlife sanc^JaIy, nature res^e 
or private reserve with a furdier 20% being held as Forest reserve. The latter is inoieasingly 
under threat fiom agricultural incursion and periodic excision by govan£se!!t. Finally, around 
1% of Beli2ae &Ils within Indian reservations. Use-rights to these lands have traditionally been 
allocated by Maya community groups. 



Table 3: Land area by tenure class 


status 


Area(Ha) Percent (%) 


Forest reserve 


431133 


20.46 


Protected areas 


312063 


I4.1S 


National land 


479478 


21.75 


bdiaa reserve 


26114 


1.18 


Private land 


89Q378 


40.39 


Total 


21S9166 


97.93 



Copv'ighted material 
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The difGcuIties traditionally associated witli communal land tenure combine with physical and 
cultural £actors to ensure that a comparatively high proportion of these lands are devoted to 
milpa. 

5.0 Data 

5.1 Data framework 

We estimate the model of equation (4) using data on a sample of points of land. This information 
is derived fix)m a GIS (Geographic Information System) database: that is, digitally coded vaaips 
representing a variety of land characteristics of interest. Sampling was performed using a 3-km 
rectangular grid (see Map 3), yielding 2401 mainland points with data on land cover. Of tiiese, 
46 water and uifoan points were excluded fiom the sample. 

To visualize the data extraction process, imagine stacking the data layers (mf^s) of interest. A 
pin pierces the stack at each sample point, and &e mapped information for the point — slope, 
distance to road, soil quaUly — is recorded and collated. (See Figure 2) In practice, we were 
concerned that the data layers might be slightiy misregistered. To minimize problems of 
misregistration, we pierced tiie data layers with a virtual apple corer rather than a virtual pin. 



Figure 2 - The Data Framework 



Point(Ay)wlth 
shared location 



GIS Data 
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That is, "we abated a 100-meter diameter buffer (circle) around each sample point on each data 
layer, and coded the characteristic which dominated the buffer . 

The data layers are described below. All were made available by the Land Information Centre 
(LIC) of Belize's MinisHy of Natural Resources. 

S2 Land use data 

Land use or land cover constitutes die model's dependent variable. The data are derived from a 
recent study (LIC, 1994) by the Lands and Surveys Department of the Ministry of Natural 
Resources with assistance £oom FAO. These land use maps, with a scale of 1 :S0,000, were 
produced tilirough interpretation of scaled false color prints of SPOT multispectral data. 
Interpretation involved considerable ground-tnithing. In order to obtain cloud-free imagery, it 
was necessary to use SPOT data ranging from 1989 to 1992, so the resulting map does not carry 
a single date. 

The map describes a total of 31 categories of land use/land cover. As in many studies of tins 
nature, difficulties occurred in defining distinct classes (particularly in the savanna/thicket/shrub 
areas) and figures given to rqnesent broadleaf forest cover are known to include a considerable 
amount of secondary regrowth (7 years and older). Often two agricultural uses would be 
intemnin^ed or be otherwise hard to distinguish. For instance, mechaiuzed amiual crops and 
pasture tended to be associated. For the pii;poses of this study, only a very aggregate 
classification was needed. Milpa and non-mechanized annual cultivation were classed as "semi- 
subsistence"; all other agriculture, together witii a small amount of land cleared but not undor 
cultivation, was classed as "commercial"; and all other areas, including secondary regrowth (but 
excluding water bodies and settlements) were classed as "natural vegetation"^ . 



8 This piecantioii was probably uimecessaiy. In tbc case of land cover, for instance. 85% of the buffers were 
homogeDous. and a single land use accounted for more dian 70% of the area in 95% of all buffers. 

9 Almost aU the "natural* v^etation in Bdize has been modified by fanmanact^ Ibe category names are for 
convenience only. 
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S3 Land systems data 

The land systems data describes the soil's physical and chemical characteristics. These data are 
taken £rom a series of Land Resource Assessments (King et ai. 1986,1989,1992), which were 
designed to yield planning information on the land's suitability for alternative crops. 

The Land Resource Assessments were based on a combination of aerial photography and field 
surveys. The methodology involves segmenting the landscape into microregions, based on 
agricultural potential as predicted from topography, soils, and vegetation. The land systems map 
is segmented into about 10,000 of these microregions, which fall into 350 distinct classifications, 
called land subunits'^. Eadi subunit is characterized by a set of physical and chemical 
descriptors, vMch in tum are used to assess suitability for ^h of 19 agricultural land uses. It is 
important to stress that nutrient values are derived fi^om field sites which are ot in agriciiltural 
use. These values are then imputed to occurrences of the same land subunit which are under 
cultivation. 

For the purposes of this analysis, it was sometimes necessaiy to assign point values to a variable 
(such as pH) whoe the land systems classification assigned a range. In other cases ~ especially 
for categorical variables ~ die values were aggregated so as to constitute binary dummies. 

Land systems variables used in this analysis are as follows. (Original coding and other 
information are in King et al, 1993, pp 110-1 17). 



10 These are similar to die land facets commoidy associated with ^ land systems methodology. 
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SLP25PL —Dummy variable, slopes 25 degrees or higher. (From a 6-category slope descriptor.) 

FLOODHAZ ~ Dummy variable for flood hazard. Floodhaz=0 for areas where floods occur less frequently than 
every 20 years, or for high floodplain bench backland. Floodhaz=l for all other areas. (Original variable: 16- 
category flooding risk.) 

SANDY— Dummy for sandy soil. (Reclassified from base variable: sandy if moisture availability > 2) 

STONY— Dummy for stony soil. (Reclassified from base variable, workability; stony if woikability = "stony or 
compact", "very stony", "25 to SO cm to bedrock", "very stony and very shallow", "stony, very shallow and 
imperfect drainage; or very compact") 

WORKABLE — Dummy for soil which is "easily woikable" (Reclassified from base variable, workability. Omitted 
categories - neither workable nor stoity, are: "rather compact", "imperfect drainage", "poor drainage", "variably 
shaUow".) 

WETNESS — An eight-point ordinal scale for drainage "as indicated by how the effect of poor drainage affects die 
growing of cacao, citrus, and pasture", ranging ftam 0- 'well-drained" to 7="permanently wet". 

WWJiJTR — Dummy for total nitrogen less dian or equal to 0.2% . 

AVAIPHOS- Available phosphorus in ppm. Bray for acid soils, Olsen for alkaline. Recoded at midpoints of 8 
category scale. 

PHLEVEL — Soil pH, recoded at midpoints of 1 0 category scale. 



SA Land tenure data 

The land tenure data is the most cuirent available. The foUovdng categories are distinguished: 
national parks, private reserves, forest reserves, national land, private land, Indian reservations. 
Areas designated national land include lands leased out Some private lands may be 
misclassified as national land. 



Copv'ighled inateriai 



24 



A Spatial Model of Land Use in Belize 



5.5 Distances to roads and market locations 

Road network data is fiom a 1:50,000 topographic map series based on I980's data, updated in 
certain areas to 1993" . For each sample point, distance to road is the straight-line distance to 
the nearest point on the road network. That "trailhead" point is then used as the basis for 
calculating the on-road travel time to the nearest town. The on-road travel time from the 
"trailhead" to each town is based on the optimal (time-minimizing) route, as determined by the 
NETWORK module of ARC/INFO, the GIS software used in this project In computing travel 
times, estimated travel speeds'^ were assigned to each class of road, as follows: 

Class Description Est, maaanmni speed (mph) 

1 Major roads (sur&ced) 65 

2 Secondary roads (surfaced/unsuifeced) 45 

3 Major tracks (unsurfaced -passable all year 35 

4 Minor tracks (unsur&ced- passable in the dry season) 20 



Towns used for distance calculation are Belize City, Belmopan, Corozal, Dangriga, Orange 
Walk, Punta Gorda, and San Ignacio. These comprise the national capital and district 
headquarters, and are the only sizeable setdemeots in the country''^ . 



11 It would have been preferable to recreate the road network as of 1989, the earliest date foe the land use data. 
This will he attenq>ted in future work, using satellite imagery. 

12 For ±e purposes of die model it is sufGcient that the inteiclass ratios of diese q>eeds he appr ox imately 
correct. In future waric we will attempt to incorporate engineering estimates or observed trafBc data. 

13 la addidon, for each sample point we calculated die <m-road distance to Pomona (site of the citms-processiag 
plants) and to the nearest sugar-processing plant (die closer of Libertad and Otangt Walk), but diese were not 
used in die final model. 
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6.0 Results of estimation 

6.1 Coefficient estimates 

The multinomial logit model was estimated for a diree category classification of land use: 
natural vegetation; semi-subsistence farming (including milpa and nonmechanized annual 
crops); and commercial farming (primarily sugarcane, mechanized annual crops and pasture, 
with some citrus, bananas and minor agricultural land uses). National parks and private 
reserves were excluded from the sample, because only one of the 322 sanqile points in these 
classes was used for agriculture'* . 

The lesults are shown in Table 4 with variable means reported in Table 5. The overall 
estimate is significant at the p=.0001 level. The coefficients for natural v^etation, the 
comparison class, are nonnalized at 0. 

Both agricultural land uses imply strong declines in rent with increasing distance from 
road. The decline is much sharper for commercial &iming, however. Similarly, on-road 
distance to market depresses rent more rapidly for conunercial fuming than for semi- 
subsistence farming (the small absolute difference in the coefficioits on distance squared 
results in a substantial difference in impact.) 



14 In effect, we assign an infinite negative coefficient to indicators for these classes, which is, fbr all ptactical 
purposes, die tnaximnm likelihood estimate. 
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Table 4 - Multinomial logit estimate of land use 

Niunber of obs b 2021 
Chi2(40) t= 671.58 

Prob > chi2 = 0.0000 
Pseudo R2 = 0.3932 







Coef . 


Std. Err. 


z 


P>| 




SBKI-SIIBSISTBNCB FARMING (MIX.PA AMD 


OTHER NONMBCBANZZBD) 








Slp25pl 


DllMinr, SLOPE > 25 


-.7597318 


.4892326 


-1 .553 


0. 


120 


floodhaz 


DUMMY, FLOODHAZARD 


.5973458 


.4083262 


1 .463 


0. 


143 


Seoidy 


DUMMY 


.9B7341B 


.5682349 


1 .738 


0. 


082 


stony 


DUMMY 


-.6283948 


1.183857 


-0 .531 


0. 


596 


workable 


DUMMY 


- .7156817 


.5879095 


-1 .217 


0. 


223 


wetness 


8 POINT ORDINAL SCALE 


-.2856389 


.1024821 


-2 .787 


0. 


005 


natlland 


DUMMY, NATIONAL LAND 


.7025545 


.3448903 


2 .037 


0. 


042 


lowjaitr 


DUMMY, NITROGEN < 0.2% 


-1.433593 


.4144497 


-3 .459 


0. 


000 


avalphos 


AVAILABLE PHOSPHORUS, PPM 


.0425456 


.0338525 


1.257 


0. 


209 


phlevel 


PH 


5 . 062912 


2.277653 


2.223 


0. 


026 


pblevel2 


PH SQUARED 


-.4143897 


.1987818 


-2.085 


0. 


037 


d_to_rd 


KM TO ROAD 


-.4435889 


.186327 


-2.381 


0. 


017 


ds_tord 


KM TO ROAD, SQUARED 


.0067064 


. 0210698 


0 .318 


0. 


750 


d3_town 


ON ROAD TIME TO NRST TOHN 


-.062215 


.0230098 


-2,704 


0 . 


007 


dSstovm 


ON ROAD TIME, SQUARED 


.000535 


.0002576 


2 .077 


0 . 


038 


X 


X-COORDINATE 


-.9751901 


.3428826 


-2 .844 


0. 


004 


y 


Y-COORDINATE 


.5322553 


.3089695 


1.723 


0 . 


085 


xy 




.6364214 


.2440257 


2 .608 


0 . 


009 


X 2 




-.3403938 


.31772 


-1 .071 


0 . 


284 






-.2012505 


.0975561 


-2.063 


0. 


039 


_C011B 




-341.2074 


256.0914 


-1 .332 


0 . 


183 


CQMHBRCZAL PAKMXNO (Sugar, BMchanlzed food crppa 


V pasture. 


citrus) 






Sip25pl 


DUMMY. SLOPE > 25 


-.9772169 


.5160684 


-1 .894 


0 . 


058 


floodhaz 


DUMMY, FLOODHAZARD 


.3284057 


.3194296 


1 .028 


0. 


304 


sandy 


DUMMY 


.8330149 


.3626363 


2 .297 


0 . 


022 


stony 


DUMMY 


-.1607761 


.4406501 


-0 .365 


0. 


715 


workable 


DUMMY 


-.1541151 


.3345063 


-0 .461 


0 . 


645 


wetness 


8 POINT ORDINAL SCALE 


-.1405025 


.0S54008 


-2 .148 


0 . 


032 


natlland 


DUMMY, NATIOIAL LAND 


-.1526676 


.2644397 


-0 .577 


0 . 


564 


lowjtiitr 


DUMMY, NITROGEN < 0,2% 


-.951181 


.2620511 


-3 .630 


0. 


000 


avaiphos 


AVAILABLE PHOSPHORUS, PPM 


.0603394 


.0224646 


2.686 


0. 


007 


phlevel 


PH 


.2163639 


1.302868 


0.166 


0. 


868 


phlevel2 


PH SQUARED 


.0230276 


. 1111243 


0.207 


0. 


836 


d_to_rd 


KM TO ROAD 


-.9344324 


.1634352 


-5.717 


0. 


000 


ds_tord 


KM TO ROAD, SQUARED 


.0480959 


.0188488 


2.552 


0. 


Oil 


d3_town 


ON ROAD TIME TO NRST TOWN 


-.0644578 


.0210767 


-3.058 


0. 


002 


dBstown 


ON ROAD TIME, SQUARED 


.0001011 


.000276 


0.366 


0. 


714 


X 


X COORDINATE 


1.414435 


.2211969 


6.394 


0. 


000 


y 


Y COORDINATE 


-.4319792 


.1107751 


-3.900 


0. 


000 






-.7925672 


.1037173 


-7.642 


0. 


000 


X 2 




.1247388 


.1833633 


0.680 


0. 


496 


Y_2 




.1822047 


.0325773 


5.593 


0. 


000 


_cons 




176.7267 


106.2071 


1.664 


0. 


096 



Note: Coefficients for contparison groiq) (nonagricultural vegetation) set to 
zero. Excluded from analysis: national parks, private reserves, towns, water 
bodies 



Multinomial ragrasBion 

Log Likelihood =^ -518.19933 
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Table 5- Variable means 



Variable 


Obs 


Mean 


Std. Dev. 


ni.n 


Max 


id 


2021 


11223 . 94 


711 . 5796 


10007 


12513 


vetoes s 


2021 


2 . 190995 


o ji o c o A a 

2 .425348 


0 


7 


floodhaz 


2021 


.4067293 


.4913451 


0 


1 


sandy 


2021 


.5299357 


m n oo O 

.4992266 


0 


1 


workable 


2021 


.2424542 


. 4286737 


0 


1 


stony 


2021 


. 0821376 


O O ^ f it o^ 

.2746421 


A 

0 


1 


low_nitr 


2021 


. 5893122 


. 4920804 


A 

0 


1 


steep 


2021 


.1459673 


O M T f O A 

. 3531609 


A 

0 


X 


stpkarst 


2021 


. U /J.2317 




A 

u 


X 


X 


2021 


320 .4o22 


O O A it O 

29 • 042 


2fi2 . 3 


Jo2 . 3 


y 


2021 


1905 . 332 


OO 1 & Of 

70 .2906 


1757 . 5 


A A Jl *> C 

2042 . 3 




2021 


. 0554181 


OO OOC*I 

. 228o512 


A 
U 


X 


d_to_ra 


2021 


3 . 086103 


3 .637712 


.0009148 


18 . 05557 


xy 


2021 


611 . bo!^? 


69 . 03314 


A C^ '3 A ^ O 

461 . 3438 


lie 








Xo . 3 /3j I 




X4D - JUO.£ 




2021 


3o35 . 989 


o^T o 1 on 

267 . 3133 


3088 . 806 


jl n T*! O Af 

4X7X . oUb 


03 covm 




O ^ Tl AO C 

J3 . / J.Uob 


. If J.41 / 




XX3 . /X 


d3_suga.r 


2021 


102 . 5999 


^o Q>f n OO 

62 . 94123 


o o 
• 28 


235 . 82 


oa cown 


2U2X 


^cf . o J9d1 




X • j9o44£b 


bo . ays / X 


da_su.gar 


2021 


102 . 9492 


64.44146 


ji >■ o n 1 no 

.4420102 


ACA CO OO 

252 . 5397 


avaipboB 


2U21 


/ . 32X324 


4 . o Jc5Uo / 




23 


nacj-xana 


2u21 




. 4U02U2X 


A 

u 


X 


ind_x'es 


2U21 


. 0J.J.337 / 


. H40J0 


A 
U 


X 


XXO. t« V CM w 


5 021 

A WAX 






0 




d3stown 


2021 


1569.437 


1871.254 


.16 


12929.96 


d3ssugar 


2021 


14486.38 


14496.31 


.0784 


55611.07 


ds_tord 


2021 


22.75043 


48.17415 


8.37e-07 


326.0036 


slp25pl 


2021 


.2172192 


.4124551 


0 


1 


phlevel 


2021 


5.715064 


.9828046 


3.5 


8.3 


plilevel2 


2021 


33.62739 


11.2711 


12.25 


68.89 


pnatvg 


2021 


.8852053 


.1791577 


.1151127 


.9999737 


pmilpa 


2021 


.0301831 


.0775559 


3.12e-08 


.7695376 


pcomfarm 


2021 


.0846116 


. 1664975 


1.58e-06 


.8846226 



These tesults ate plausible. Commercial burners are likely to be higUy sensitive to road and 
market access. This vnil be particularly true for sugar &im«s, viiere typical transport costs 
range &om$10-$18 per ton, against mill prices of $55. (King etal., 1993) Milpa and other small 
&nners may market only a fraction of their produce, making them less sensitive to distance to 
market and to road. For ccops tiiey do market, transport costs may be less importanL Typical 
transport costs for dryland maize in Toledo district, for instance, are just 5.6% of mill price'^ 



15 Of coarse, tiic transport cost rado nu^ increase dramatically with off-road distance. 
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(King et al., 1993). Nonetlieless, %mi-subsistence fanners may value access to schools, clinics, 
and off-fann employment sites. 

The predicted iiopact of soil chemistry also appears to be plausible. Nitrogen has a strong, 
hi^y significant effect It is relatively more important for semi-subsistence fumers, either 
because of differing crop mixes or because of credit constraints in purchasing fertiliz^'^ . Soil 
acidity (as measured by pH) has a strong, highly significant negative effect on bid rent for semi- 
subsistence Arming. 

The log odds in favor of semi-subsistence farming increase by 2.8 as pH increases from its 
mimmiim of 3.5 to an Optimum at 6.1; above that point, more alkaline soils are predicted to be 
less desirable. Soil acidity also depresses rent for commercial &iming (the two terms are jointly 
significant at the .05 level), but there is no penally for alkaline soils. Finally, the point estimates 
for the effect of available phosphorus are positive — highly significant for commercial Arming, 
not significant at conventional levels for semi-subsistence fiuming. In general, these results are 
consistent with anthropological findings that traditional &nners use a variety of pedological and 
botanical cues to assess soil quality vdth considerable accuracy (see Carter 1969 on Maya 
fanners in Guatemala; Wilken 1987 on Mexican and Guatemalan farmers; and Moran 1993 on 
Braalian &rmers). 

Turning to o&er land characteristics the effect of high slopes on rent is negative for both land 
uses, as expected. The coefS^cient for semi-subsistence farming is just shy of significance at the 
.1 level and is smaller in absolute value than that for conunercial farming. While far from 
conclusive, this is consistent wi& hilly land being relatively better suited for nonmechanized 
than mechanized farming. Wetness has a substantial and highly significant negative effect on 
TexA for both land use categories. The effect is however twice as strong for semi-subsistence 
Arming, perhaps because wetter soils are relatively harder for nonmechanized &rmeis to work. 



16 AltenKttively , this nu^ a clue that soil nutrient deficiencies need to Iw inimcted widi distance in ttie 
equatiQI^ it is more expensive to remediaiB poor soils when th^ are fiuiber fmm. the road. 
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Indicators for sandiness and flood hazard bad estimated positive effects; sandiness was 
statistically significant. These variables may reflect tiie desirability of land near rivers. An 
indicator for stony soil had a negative coefficient, but, surprisingly, was not statistically 
significant The only tenure variable included was national land, which had a strong positive 
effect on tiie value to semi-subsistence &rming, as expected. 

Finally, the spatial trends play a major role in fitting the model. For semi-subsistence Arming, 
the standard deviation of the implied spatial effect is 2.1, with a range of 7; for commercial 
timing the standard de^dation is 1.28, witii a range of 3.2. The spatial effect for semi- 
subsiscoice &rming has its peak in the Mayan Indian areas of Toledo district The corresponding 
peak for commercial Arming is near the sugar processing plants in the north. Linear and square 
terms in distance to the nearest sugar processing plant were jointly insignificant when added to 
the model because of high multicollinearity with the spatial trend variables. 

To check the robustness of the multinomial logit model, binary logit and probit models were run. 
In these models, the two classes of agricultural use woe combined, and additional e^lanatoiy 

1 ^ 

variables were used . The results for both specifications are qualitatively very similar to die 
multinomial logit results. Table 6 shows the binary logit results. In the binary model, it is 
possible to disaggregate high-slope land into non-kaist and karst; the latter much more strongly 
discourages cultivation. It is also possible to add a dummy variable indicating high altitude 
(MOO meters). This too strongly discourages cultivation. The two tenure variables — indicators 
for Indian reservation and for national lands — have positive coefficients but are not statistically 
significant The modest and insignificant net impact of these variables reflects strong 1wt 
opposing impacts on subsistence versus commercial cultivation. 



17 Dummy variables which perfecdy predict one but sot both land uses (e.g., there is some subsistence 
agriculture but no commercial agricultuie in Indian leservadons) present estunation problems for standaid 
software in the multinomial case, but are easily accomodated here. 
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Table 6: Binary logit estimates 

Dependent variable: Agricultural use (n= 2021) 



Parameter 


QtnnHfifvl 

OuUlwUU 


WalU 


Pr> 




Variable 


DF 










INTERCPT 


1 




140 1 






HI 


1 




1 0789 


3 7512 


0 0528 


STEEP 






0.41 56 


1 0756 


0J2997 


STPKARST 


J 




n S44n 


4 0A5fl 


0 0259 


FLOODHAZ 


I 


0.4384 


0.2519 


3.0299 


0.0817 


SANDY 


1 


0.6388 


0.2861 


4.9847 


0.0256 


STONY 


1 


-07S60 


0 4145 




0 0681 


WORKABLE 


1 




0.2787 


2.4660 


0.1163 


WETNESS 


1 


-0.2235 


0.0542 


17.0360 


0 0001 


LOW NTTR 


I 


-0.8151 


0.2165 


14.1778 


0 0002 


AVAIPHOS 


1 


00620 


0.0179 


12 0346 


0.0005 


PHLEVEL 


1 


-0.0551 


1.1239 


0.0024 


0.9609 


PHLEVEL2 


1 


002S8 


0.0962 


0 0720 


0.7885 


IND RES 




03969 


0 5283 


1 2767 


0 2585 


NATLLAND 


J 


02007 


0.2112 


0.9031 


0.3419 


D_TO_RD 




-0.7099 


0.1181 


36!l556 


0.0001 


DS TORD 




0.0289 


0.0122 


5.6291 


0.0177 


D3 TOWN 




-0.0503 


0.0155 


10.5220 


0.0012 


D3ST0WN 




0.000198 


0.000195 


1.0366 


0.3086 


D3 SUGAR 




-0.0341 


0.0116 


8.6042 


0.0034 


D3SSUGAR 




0.000104 


0.000049 


4.5577 


0.0328 


X 




0.5391 


0.1435 


14.1144 


0.0002 


Y 




-0.4027 


0.1555 


6.7100 


0.0096 


XY 




-0.4173 


0.0845 


24.4117 


O.0OO1 


X 2 




03660 


0.1660 


4.8613 


O.0275 


Y 2 




0.1386 


0.0421 


10.8563 


0.0010 



6.2 Predictive ability of the model 
6.1 Spatial patterns 

Maps 4 and 5 show the predicted probabilities of agricultural use over 1lie landscape, for semi- 
subsistence and commercial agriculture. The probability peaks correspond reasonably well with 
the actual land use patterns shown in Map 2. The predicted commercial agricultural areas 
coincide wilfa the sugarcane region of northern Belize, the mechanized &rms of the nor&\\^ 
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and center-west (near San Ignacio), and the citrus-growing areas near Dangriga. The predicted 
semi-subsistence areas are centered on the actual subsistence farms of Toledo district, in the 
south. 

Note that the model tends not to predict particular points as having very high probabilities of 
agricultural use; rather, it predicts entire areas as being predisposed, with each individual point 
having a modest probability of cultivation. Is this fuzziness a vice or a virtue? It is a vice if it 
reflects omission of important information which differentiates the agricultural suitability of 
neighboring points. It is a virtue if neighboring points are in fad very similar in agricultural 
suitability. In that case, it is truly a matter of chance which points are currently under 
cultivatioiu'^ This is particularly true for shifting cultivation, where fields are rotated between 
forest Mlow and crops. 

Classification accuracy 

Because we argue that even low predicted probabilities of agricultural use convey information, 
we use a 20% probability as the prediction threshold (Le., designate a point as predicted to be in 
agricultural use if tibe predicted probability is greater than 20%). Using tins criterion, tiie model's 
accuracy is 76% for actual agricultural points, and 89% for actual natural vegetation points 
(Table? 



18 The points are not all simultaneously under cultivaticm because that would diange the macro spatial price 
equilibrium. 

19 These and subsequent pre^ctive assessments exclude natural paiks and private reserves. If we included these 
points and followed die rule that all such points at? :;iiedicted to be under natural vegetation, the overall model's 
predictive performance would be mudi improved. 

20 Incontrast. xise of a conventional S0% threshold dnunaticaSy reduces accuracy £orcdtivaledp(»nts^ 
and boosts accuracy for noncultivated only modesdy, to 98%. Indeed one can trivially achieve an overall 
accuracy of 88.5% sinq>ly by predicting all points to be noncultivated. This illustrates the familiar statistical 
trade-off between type I and type 11 errors, and die artntrariness of die 50% direshold. 
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Table 7 - Predicted vs. actual agricultural land use 





Prob(agric.) <2 


Prob(agric.) >.2 


TOTAL 


Actual natural 
vegetation 


1591 


198 


1789 


Actual i^cultural 
use 


56 


176 


232 


TOTAL 


1647 


374 


2021 



The next table (Table 8) cross-tabulates actual land use against predicted use for semi- 
subsistence agriculture, again using a 20% probability threshold. What is particularly interesting 
is that of the 38 natural vegetation points incorrectly predicted to be in agriculture, 31% were 
secondary regrowth or thicket, indicating recent cultivation. By contrast, only 6.5% of the 
correctly predicted natural vegetation points were secondary regrowth or thicket This further 
supports the assertion diat low predicted probabilities of cultivation are usefid for prediction. 
Note that the prediction criterion distinguishes sharply between conunercial aiMl semi-subsistence 
agriculture. 



Table 8 - Actual land use vs. predicted probability of semi-subsistence agriculture 





Prob(sani-sub) < .2 


Prob(seinj-sub) > 2 


TOTAL 


Actual nstutal vegetation 


1751 


38 


1789 


Actual semi-subsistence 


32 


29 


61 


Actual sugarcane 


64 


0 


64 


Actual odier commeicial 
agriculture 


103 


4 


107 


TOTAL 


1950 


71 


2021 



Finally, the following table (Table 9) compares predictions for commercial cultivation with 
actual land use again using a 20% threshold. Fredictioiis are highly accurate for points under 
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sugarcane, less so for other commercial agriculture. Semi-subsistence and commercial 
agriculture are sharply distinguished. In fact, only one sample point had predicted use 
probabilities of more than 20% for both semi-subsistence and commercial farming. 



Table 9 - Actual land use vs. predicted probability of commercial agriculture 





Prob(commercial 
agriculture) < .2 


Prob(cominacial 
agriculture) >.2 


TOTAL 1 


1 Actual natural vegetation 


1661 


12S 


1789 


1 Actual setni-subsistence 


54 


7 


61 


Actual sugarcane 


3 


61 


64 


Actual other commercial 
agriculture 


43 


64 


107 


TOTAL 


1761 


260 


2021 



7.0 Discussion 

7.1 Implications for road-building 

The results show that distance to road, on-road travel time to market, and agricultural suitability 
have a strong impact on land use and deforestation. 'Hiese results are illustrated in Figures 3 and 
4, vMch show predicted land use as a function of the three variables. Figure 3 describes land 
use on "agriculturally suitable" soils, here defined as not high slope, not low nitrogen, and sandy, 
with other soil \^riables set at the conditional meaiL Spatial trend effects for both agricultural 
land uses are set at sample means. The three panels p ^sent predictions across two dimensions of 
distance to market: the distance to the road Qiorizontal axis of each panel) and the subsequent on- 
road distance to the nearest town (differentiated by paiiel). The top panel describes points -wlaciL 
are close to town: the on-road travel time is just one minute. For points which are on the road 
(distance to road=0), Hxic is a 62.4% probalnUty of agricultural use, consisting of a 49.2% 



Agriculturally suitable land 

On-road time to town = 1 minute 



Figure 3 




On-road time to town = 20 minutes 




0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0 



4 5 6 
1^ to ro^ 
On-road time to town = 1 00 minutes 



10 



ISemi-subsistence 
ICommerdal 




4 5 6 
Km to mad 



8 



10 



Copv'ighlec) inaterial 



A Spatial Model of Land Use in Belize 



37 



probability of commercial agricultural use, and a 13.2% probability of semi-subsistence 
cultivation (in this, likely some form of truck fimning.) As distance to road increases, the 
cultivation probability decreases, witii commercial farming failing off more rapidly than 
subsistence farrtaag. At 10 km from town, tiiere is a 97.8% probability that the land is still 
under natural vegetation. The middle panel repeats the predictions for points sharing an on-road 
travel time of 20 minutes. The agricultural use probabilitiea are substantially reduced. The 
bottom panel describes points with an on-road travel time of 100 minutes. At this distance, the 
probability of commercial use is nearly zero; the land is more valuable for semi-subsistence 
fiirmers. Even for these fanners, the probability of cultivation fidls off rapidly with distance 
from the road. 

Figure 4 repeats the set of predictions for agriculturally marginal land, defined as being high 
slope and low nitrogen, witii other variables set at their conditional means. The contrast with 
agriculturally suitable land is striking. Even close to town (top panel), the probability of 
agricultural use is just 15.2% at the road, dropping below 1% at S km from tiie road. At just 50 
minutes (on-road) to town, the probability of agricultural use is only 1.2% at the road, dropping 
to 0.6% 1 km off the road. The dominance of commercial over semi-subsistence framing results 
from the latter's much stronger disutility for low-nitrogen soils. Clearly, though, these results are 
shaped by Belize's low population density. 

With these estimates, we can begin at least indicatively to assess the costs and benefits of roads 
intensification vs extensification strategies. The top panel of Table 3 suggests high economic 
returns to inopeasing tiie density of roads in frtvorable areas near market c«itets. Building a 
feeder road to points which were formerly 2 km fix)m die road ^stem, for instance, boosts 
cultivation probability from 28.2% to 62.4%. The impact on labor donand and by extoision on 
poverty dq)ends on the particular crops and technology chosen. We have not modeled 
agricultural intensification here exphdAy, but theory (see Appendix A) and expenence suggests 
more land and labor-intensive agriculture nearer to town. The environmental impacts of such an 
intensification strategy will vary with tiie local ckcumstances, but in geooal we 
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would expect areas closer to town to have lower value for biodiversity conservation, simply 
because they are more ^sturbed and fragmented. 

A contrasting, extensification strategy would push roads out to more remote areas. Newly 
opened areas vdth good soils would experience some forest conversion to semi-subsistence 
cultivation — in this case, almost certainly milpa l^ucming. Returns to labor and land would be 
relatively low. The environmental consequences might be substantial. Crop rotation woidd 
affect an area several times larger dian that indicated by the current cultivation probabilities. 
Habitat fragmentation might result If the roads were extended into remote areas with poor soils, 
there would be almost no economic gains to counterbalance the ex^pexxses of road-building. 
While there would also be litde impact on forest cover, road access would e^ose the forest to 
various forms of degradation, such as overextraction of mahogany or poaching of birds. 

We stress that this is only an uidicative analysis. A more thorough analysis would involve 
calculation of the impacts of particular road-siting alternative. It would also allow for general 
equilibrium efifects: a substantial increase in cropped area would boost wages and reduce die 
price of domestically-consumed agricultural products, changing the coefficients embodied in the 
model. Finally, it would take into account distributional effects (across income groups or 
regions) of altered cultivation patterns. 

7.2 Effectiveness of habitat protection 

For conservation planning, it is useful to be able to a^ess the effectiveness of habitat protection. 
In the case of Belize, we observe very little cultivation in national parks and private reserves 
(about 0.3%). Is this because these areas are effectively policed, or is it because they are remote 
or otherwise unattractive for cultivation? 

To address tbis question, we use the estimated coefGcients to predict the extent of cultivation 
tibese areas would experience if they w»e not jirotected. Taking die expectation of predicted 
probabilities, 1.3% of the area would be predicted to be under current semi-subsistence 
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cultivation, and 2.7% under commercial cultivation. Since the area affected by shifting 
cultivation could be five to ten times larger than that under cultivation at a given time, we 
conclude that habitat protection has been effective in Belize. 

7.3 Effect of logging 

A striking finding is the predicted very low probability of agricultural use in areas which have 
poor soils and are far from madcets — even if they are close to roads. In fact it is easy to identify 
areas which are close to main roads but are not used for agriculture (see Map 6). This suggests 
that the hypothesis (ofien applied to Asia) that logging causes damage primarily by inducing 
follow-on migration does not necessarily apply to low-population density, remote, areas. In light 
of the results presented here, this suggests a more detailed examination, ideally using GIS data, 
of the dynamics of logging and deforestation in areas fai &om markets. 

8.0 Summary and next steps 

We have presented a static equilibrium model which relates observed land use to the relative 
returns to different alternatives. Relative retums at each point in the landscape are determined by 
road access, distance to market, and the inheroit productivity of the land. The model was 
strongly supported by data fiom Belize. Market access and distance to road strongly affect the 
probability of cultivation; as hypothesized, the effect is stronger on commercial than on 
noncommercial agriculture. High slopes, poor drainage, and low soil fertility discourage both 
commercial and semi-subsistence agriculture. Semi-subsistence {^culture is found to be 
particularly sensitive to soil acidity and lack of nitrogen. This confirms anthropological studies 
which find that subsistence burners are shrewd judges of soils. Taken together, these results 
suggest that intensification of tiie road network around market areas offers higher economic 
retums and lower environmental impacts than extensification of the road network into new areas. 
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The model has a Avide variety of shortcomings which we hope to conect in future work. First, it 
treats load luting as exogau>iis. Although data on agricultural suitability helps control for any 
resultant biases, we have suggested instrumental variables strategies to address this problem. 
Second, the assumptions tmderlying the multinomial logit specification are very strong; they will 
be relaxed via die use of multinomial probit or nested logit specifications. Third, the spatial 
trend variables currently explain much of the spatial variation in land use - \^ch is to say, the 
underlying causes have not been identified. We hope that incoipoiation of information on ibc 
spatial distribution of population will improve die model's explanatory power. 

Two immediate extenaons are planned. First, we will examine the alnlily of the model to predict 
land use changes over die period 1985-1994 in the Toledo District, die site of rapid recent 
immigration. Second, we will assess the potential impact of road network extensions on the 
environment, using data on die spadal distribution of species. Ova die longer run, we would 
like to e}q)licidy introduce spatial price formation into die model, so as to be able to eiq^lore die 
implications of changes in agricultural markets and polides. 
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Appendix A: Model derivation 

Note: The equation numbers in this derivation parallel those in the non-technical presentation in 
the main body of the text 

Let be the price of the ou^ut of use k at point i^ Cjj^ be a vector of prices of inputs to k at that 
point, Xj|]r be the optimal quantities of inputs for k per unit of land, Qj^ be the potential output 
of k at the pomt, and Uj^ a random disturbance. The potential rent associated with devoting ihc 
point to use k is: 

P,C, and Q are endogenous, and in the current example unobserved. We derive a reduced from 
by specifying observable determinants of these variables. In the classic von Thilnen 
interpretation, P and C are both strongly related to distance to market (or port or processing plant, 
dependmg on the commodity in question), reflecting transport costs. While a simple von Thflnen 
model uses a linear function of distance, truncating vAiea P=0, here it is more convenient to 
specify exponential functions: 

(2) Pik=exp[yok^yjifi,^y2i;rj 
Ca=e:qj[8ok+ Bj^Pi+SfTJ 

where Dj is the distance from point i to the road, T| is on-road travel time to the nearest market, 
and the parameters are commodity-specific. (For a closer fit to a truncated linear function, 
polynomials in D and T can be added.) We hypothesize HaaX output price decreases with both 
dimensions of distance (y^>0, y>0)j and that input costs ina:ease with distance from market 
(5>0.8 >0). 



The production fimction, here e>qpressed as output pex unit of land, is assumed to be: 
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(3) Q.k=s.^r [0<M1] 

The productivity factor S can be expressed as the product of agroclimatic and other variables si, 
S2,... describing land characteristics: 

Ik ~ 'X)*«^I/ "2/ ■■• 

From (3), we can derive the demand for X: 

(.) ^= [^]'^« 

Substituting (3) and (*) into (1), suppressing the enx>T torn and the subscripts, and simplifying, 
we have: 

R=PQ-CX=FSJ^-CX=J£{PS)^^-C]= 
(♦*) C*'^P-'^(FSp]''^'~P^(l-P)/p 

Together, and (**) show lliat rent and input-intensity increase as output price increase, and 
decrease as input costs increase. 

Substituting (2) and (3a) into (**), taking logs, reintroducing the commodity subscript and the 
stochastic term, and relabeling coefiBcients: 

(4) ln«ifc==aoifc+a,ilnD/fa2tlii7;^xj3tlnsu+...+Urt = ZjAjt+ua 

The coefficients on distance are predicted to be negative; those on productivity-enhancing land 
characteristics s are predicted to be positive. The magnitude of those coefficients will however 
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di£fer maikedly fix>m crop to crop. For instance^ transport costs are relatively high for sugar, 
yielding a large (negative) coefiQcient on the distance measures. By contrast, milpa fermers viio 
produce mosdy for sub^stence may care relatively litde about distance to market 

Point i is devoted to use k if: 

(5) R^R^ 2Xij*k 

Returning now to the stochastic specification, if the disturbances u are Weibull distributed and 
uncorrelated across uses j, then tiiis is equivalent to a multmomial logit model where the 
probability that plot i is devoted to use k is: 

J 

(The coefficients of one use must be normalized to zero. Here forest is a nabiral choice for the 
default cat^ory.) 
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Appendix B: Instmmental variables approaches 

Suppose that road siting is detennined, in part, by indicators of agricultural development 
prospects. Suppose, fiirlher, Ifaat we lack information on those indicators. Then our estimates of 
the efifect of road proximity on agricultural development will be biased, since proximity and 
development are co-determined by unobserved variables. 

The instrumental variables solution is to find variables which are effective at explaining distance 
to the road or on-road travel time, but vWiich are not correlated with agricultural development 
(after controlling for observable land characteristics). Auxiliary equations are then used to 
predict distance to road and on-road travel time as a function of exogenous variables. The 
predicted distance values are used in the main land-use equation estimates. 

It is difficult to identify suitable instrumental variables. We tried an approach based on "virtual 
roads"^^ . The underlying premise is that the locations of major tov^ are exogenous, detennined 
by geography and historical accident We then argue tfiat the location of tiie towns points 
predetennines the geometry, but not the precise routing, of the main intertown road netwoiic. In 
other words, transportation demand will inevitably (aeate links between major centers, and those 
links are exogenous to the areas they traverse. 

Road network geometry is used to g«aerate instrumental variables as follows. Using GIS 
technology, we draw straightline virtual main roads between the seven market points, mimicking 
die actual primary road network. We then construct a virtual path linking any given sample 
point to the nearest location on the virtual road systeirL We can compute the mean slope and 
mean swampiness along tiiis path. The length and characteristics of the virtual path, together 
with the on-virtual-road distance to the nearest town, are used as instruments for the actual on- 
road and off-road distance between san:q)le point and nearest tovm. 



21 This strategy was wcffked out with Peter Orazem. 
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In practice, we found that ifaese instruments worked very well in explaining on-road travel time. 
In fact, straight-lme distance between the sample point and the nearest town, by itself, had very 
strong explanatory power and was fer simpler to calculate. Unfortunately, the instruments had 
litde esqjlanatory power for distance to the road. The resultant predicted values of distance to the 
road were highly collinear Avilh other determinants of land use and therefore could not be used in 
the land use equation. 

In future work, we will try an alternative approach. We hypothesize that distance to the road is a 
function of a point's accessibility, and of the agricultural suitability of nearby points. Holding 
constant the characteristics of a given point, it will tend to be fardier £rom a road if it is in a 
mountainous area or in a swamp than if it is in a fertile plaiiL Thus we can construct instruments 
by computing Ihe mean slope, elevation, swampmess, and agricultural suitability of the area 
within, say, a five kilometer radius of each point (Note that tiie land use equation must then 
include not only ciistance from the road, but also mean slope and swampiness on the shortest path 
from point to road.) In the Belize context, the local density of mahogany and cedar trees is an 
additional, potentially important determinant of road construction, though data are not readily 
available. 
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